Spin susceptibility and the 7r-excitation in underdoped cuprates. 
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The dynamical spin susceptibility x'l a t wave vector (it, tt) and the spectrum tt" of the spin-triplet 
particle-particle excitation with center of mass momentum (n, it) (7r-excitation) are considered in the 
slave-boson formulation of the t-J-model. Propagators are calculated in a diagrammatic t-matrix 
approximation in the d-wave superconducting state for a wide doping range. The resulting spectra 
x" and 7r" both show a resonance at a doping dependent energy, in qualitative agreement with recent 
numerical cluster calculations. In underdoped systems, the peak position is comparable to that found 
in neutron scattering experiments. The peak in x" as well as tt" is at low doping entirely caused by 
spin fluctuations, whereas the triplet particle-particle channel does not contribute as a collective mode. 



The spin-triplet particle-particle excitation ( l 7r- 
excitation') 

7?+ = ^(cos(fe x ) -cos(fcj,))cL feT 4 +9T 



at wave vector q — (tt, tt) has been introduced in Q as a 
possible explanation for the '41 meV resonance' observed 
in neutron scattering on cuprates in the superconducting 
state (see e.g. (|||). It has been argued that tt^ is an 
approximate collective cigenmode of the t-J or Hubbard 
model. The coupling of spin-triplet particle-particle ex- 
cited states and spin-singlet particle-hole states in the 
superconducting phase then should lead to a resonance 
in the susceptibility x"( w ) at g = (tt, tt) at the energy loq 
of this 7r-mode. 

We compare the susceptibility and the propagator of 
the 7r-excitation in a slave-boson theory for a wide range 
of hole concentrations (doping) . Both calculated spectra 
X" and 7r" show a pronounced resonance at the same en- 
ergy loo , which is roughly given by the chemical potential 
H as wo ~ 2\fi\ . The outcome is in qualitative agreement 
with the aforementioned prediction and with recent nu- 
merical j|,U and diagrammatic calculations || . How- 
ever, our interpretation differs from that originally en- 
visioned in [0 : The diagrammatic spin-triplet particle- 
particle channel does not contribute as a collective mode 
to x'i or 7T " ■ I n underdoped systems not far from the 
transition to the Neel state, the resonance is solely caused 
by the 'RPA channel', which describes spin fluctuations 
mediated through spin-singlet particle hole excitations 
of fcrmions. 

The '7r-propagator' is given as 



ir(w) = (T t t}(t)^(t')) u 



(1) 



We start from the t-J-model and consider a Gutzwiller- 
projected 7r-propagator, Eq. ^ with tt — > PqttPq . The 
calculations for tt(uj) as well as the susceptibility Xtt(u) 
are performed within the standard slave-boson scheme. 
Diagrammatic expressions are based on a self-consistent 
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FIG. 1. a: Vertex renormalized mean-held susceptibil- 
ity (part of the t-matrix approximation, see text), b: Vertex 
function for fermions in the d-wave pairing phase. The dashed 
line represents the nearest neighbor spin and density interac- 
tion ~ J of the t-J-model. 



perturbation theory with self energies taken at Hartree- 
Fock (mean-field) level. We consider the superconducting 
state at very low temperature, i.e., the d-wave pairing 
phase of fermions and fully condensed bosons. The t- 
matrix approximation for Xtt(^) and tt(lu) are indicated 
in Figs. |l| and || respectively. 

The susceptibility is given by the vertex-renormalizcd 
mean-field bubbles displayed in Fig. |l|a , which are to be 
inserted into 



(2) 



Eq. (||) represents the particle-hole RPA channel (random 
phase approximation). 

The single and double arrowed lines in Fig. |l] stand 
for the normal and pairing Green's functions of auxiliary 
fermions. The dashed line in Fig. [l]b is the t-J-modcl's 
spin and density interaction for fermions on two nearest 
neighbor lattice sites 



J [ s i s j - \nirij] 

<i,j> 



(3) 



with ni = Yu<,flf° 
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FIG. 2. t-matrix approximation for the 7r/-propagator. 
The box indicates the phase factor cos(fcir) — cos(fc M ) , k is 
the summed (loop) wave vector. The double dashed line 
stands for the effective interaction J(q,u>) (see text). Pref- 
actors count degenerate exchange parts. 

The vertex corrections entering x-n consist of the spin- 
singlet particle-hole (ph) ladder diagrams shown in Fig. 
|l|b . The double arrowed (anomalous) Green's function 
introduces the ph channel in both time directions. In 
general it also allows for a coupling of the spin-triplet 
particle-particle (pp) channel into the singlet ph correla- 
tion function Xn by transforming e.g. a spin-up fermion 
into a spin-down hole and vice versa. However, the pp 
channel would appear in x-n as a vertex-function, involv- 
ing at least one interaction vertex Eq. (||) with equal 
spin on both sites, which is zero . This reflects the 
fact that Pauli's principle blocks any exchange process 
~ J = At 2 /U for particles with equal spin in the Hub- 
bard model. Thus the pp channel contributes no spectral 
weight to x-n ■ This also holds if the pp channel is 'ar- 
tificially' switched on by replacing rurij in Eq.([| with 
griiiij and turning g = 1 — ► g = 0: Recent numerical 
cluster calculations for the t-J-model || show that x" 
is not affected by varying the coefficient of the density- 
density interaction. In the following we stick to the case 
<?=!■ 

Numerical calculations in the t-matrix approximation 
are performed with mean-field parameters set to reflect 
a fermion bandwidth of 4 J and a superconducting gap 
Aq = 40meV « 0.3 J. As has already been observed in 
earlier RPA calculations || , an instability to the Neel 
state occurs at an unphysically high hole concentration 
(doping) x c . Since the vertex corrections of the t-matrix 
approximation turn out to have no significant effect, we 
assume a further renormalization of J — » a J in Eq.(^) . 
We have chosen a — 0.5 such that x c is reduced to k, 
0.02. 

Results for x^( w ) are shown in Fig. |^(top) as con- 
tinuous curves for several hole densities in the under- 
doped regime. The dominant feature is apparently a 
sharp and strongly doping dependent resonance. Its 
position loq shifts from w at the magnetic instabil- 
ity (x = x c = 0.02) to higher energies with increased 
doping, crossing the anticipated value 40meV f» 0.3 J 
around x — x m = 0.12. ujq is for x > x c roughly given 
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FIG. 3. Top: Susceptibility x" in units of (g/i B ) 2 /(2J) 
in the superconducting phase (An = 0.3) for hole concentra- 
tions x — 0.02, 0.05, 0.1, 0.15 fom left to right. Dashed curves 
are calculated from the bubble diagrams shown in Fig. [j]a 
(multiplied by 10). Continuous curves result from the full 
t-matrix approximation, e.g., are calculated with the renor- 
malization through the spin-fluctuation (RPA) channel Eq. (^) 
taken into account. Bottom: 717-propagator in units of 1/J 
for the same set of x from left to right. Dashed lines result 
from the 1st and 2nd (bubble) diagram in Fig. ^ . Continuous 
lines include the RPA channel (3rd diagram). 

by the chemical potential fj, as w 2\fi\ . This shift of 
the resonance with hole concentration is found in neu- 
tron scattering experiments on optimal ^,|| and under- 
doped YBCO-compounds PJlO| . The 'optimal' doping 
x m « 0.12 found here also compares to experimental val- 
ues. However, the spectral weight J d?qx"(q,w)/ f d 2 q 
comes out too small with respect to the experiment ]ll[ . 

The resonance is caused by the spin-fluctuation RPA 
channel Eq.(||): For comparison, dashed curves in Fig. 
H (top) show the results for x"(w) , i.e., the renormal- 
ized bubble diagrams in Fig. [j]a . The position of the 
doping dependent peak here is bound from below by 
w 2Ao = 0.6 J , even for lowest doping. The contribu- 
tion from the ph vertex corrections is quite small, x'n 
differs only slightly from the well known mean-field sus- 
ceptibility. 

The results for the susceptibility may be compared to 
the 7r-propagator Eq. (|I]) . In slave-particle formulation, 
with bosons completely condensed at very low temper- 
ature, it reads tt(u>) = x 2 TTf(cu). The prefactor x 2 is 
the (mean-field) probability of finding two empty lat- 
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favor of the fermi-liquid state Ao = is ignored for the 
moment). The peaks in x'l an d 7r" still occur at the 
same position m 2\fi\ , but the RPA induces only a slight 
shift, besides an enhancement of x'l ■ In the normal state 
Ao = the resonance in x" vanishes, whereas the peak 
in ir" remains as a delta function , it"(uj) ~ x 2 8(oj + 2/z) , 
as has also been observed in numerical calculations on 
highly doped clusters (referenced in (l2|). In contrary 
to the low doping region, the highly overdoped regime 
» Ao is well described by mean-field theory. 
Acknowledgments: Discussions with W. Hanke and 
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FIG. 4. Set of curves for a large chemical potential /j, — — J 
in the superconducting state (see text), x'l an d 7r}' (scaled 
xl/4) are shown with the RPA-channel taken into account 
(continuous and dashed line respectively) and with the RPA 
channel omitted (dotted and dashed-dotted line resp.) . 



tice sites when adding a spin-triplet pair of particles. 
The 7r-propagator for fermions 7r/(w) appearing here is 
formally identical to Eq. (IT) , its t- matrix approxima- 
tion is displayed in Fig. E . According to the discus- 
sion given above, the triplet pp channel may contribute 
only as the mean- field bubble (1st diagram). The sin- 
glet ph channel appears in vertex renormalizations in 
the 2nd and 3rd diagram. The 3rd diagram contains 
the contribution from the RPA channel via J(q,ui)/2 — 
\J{q) + \ J {q)x{q>u)J{q) , indicated as a double dashed 
line in Fig. @. 

The resulting spectrum tt'Uuj) in the underdoped 
regime is shown in Fig. |^ (bottom) for the same set 
of parameters and hole densities x as the susceptibility. 
Continuous lines correspond to the t-matrix approxima- 
tion, dashed lines are calculated with the 3rd diagram 
in Fig. H (the coupling to the RPA channel) ignored. 
Again, the effect of the vertex corrections is negligible, 
the dashed curves differ only slightly from the mean-field 
theory (given by the 1st diagram in Fig. ||). Apparently 
7r" shows a pronounced peak which occurs at exactly the 
same position as the resonance in \'i > if the same approx- 
imation is used for both quantities. As has been pointed 
out, the spin-fluctuation (RPA) channel has to be taken 
into account in the underdoped regime < Ao , where 
the system is not far from the instability to the Neel 
state, and the RPA dominates x" . In this case the peak 
in 7r" is entirely caused by the coupling to spin fluctua- 
tions through J . Note that its spectral weight decreases 
with reduced x , and vanishes at the transition to the 
Neel state (x = x c w 0.02) . 

The picture changes in a highly overdoped situation 
3> Ao : Fig. [| shows curves for a large chemical po- 
tential ji = — J (the breakdown of superconductivity in 
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